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ON THE ANGULAR VARIATION OF SOLAR REFLECTANCE OF SNOW 


B, J. Choudluiry 
and 

A. T. C. Chang 

/ 

ABSTRACT ^ „ 

Spectral and integrated solar reflectance of non-hoinogencous snowpacks arc derived 
assuming surface reflection of direct radiation and subsurface multiple scattering. 
For surface reflection, a bidirectional reflectance distribution function derived for 
an isotropic Gaussian faceted surface is considered, and for subsurface multiple 
scattering, an approximate solution of the radiative transfer cuqation is studied. 
Solar radiation incident on the snowpack is decomposed into direct and atmo- 
spherically scattered radiation. Spectral attenuation coefficients of ozone, carbon 
dioxide, water vapor, aerosol and molecular scattering are included in the calcu- 
lation of incident solar radiation. Illustrative numerical results are given for a case 
of North American winter atmospheric condition. The calculated dependence of 
spectrally integrated directional reflectance (or albedo) on solar elevation is in 
qualitative agreement with available observations. 
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I. Introduction 


Fast growth in world population imposes higher demands on available water supplies. To meet this 
increased demand, existing water supplies will have to be managed more effectively. Runoff from 
melting snow provides greater than 65% of the total streamflow across most of the mountainous 
Western United States, In the upper midwest of United States and the Canadian high plains, know- 
ledge of snowpack characteristics is extremely valuable for flood forecasting purposes in the early 
spring. In addition, snowpacks also play an important role in climatic changes (References 1, 2, 

3 and 4). 

The variations of solar radiation absorbed and reflected by the Earth are key factors in the under- 
standing of climatic change. In a simple energy feedback mechanism of climatic change, if solar 
radiation reflected by the Earth were increased, a trend towards the extension of the winter season 
will result (Budyko-Kukla type mechanism of ice cap, Reference 5). In this context, the role 
of snowpack need not be overemphisized because it is one of the most reflective naturally occurring 
materials. An understanding of visible near-infrared solar reflectance of snow is needed for mapping 
of snowcovered areas via remote sensing methods, and estimating the solar energy reflected by the 
Earth. 

The complicated processes of reflectivity by a medium such as water, sand, vegetation or snow 
depend upon surface and subsurface radiative scattering and absorption characteristics. In general, 
the reflected radiation can be decomposed into two components: ( 1) specular and (2) diffuse. 
Radiation reflected in accordance with geometrical optics is the specular component, while radia- 
tion reflected according to the Lambertian law is the diffuse component. The relative contribution 
of these components varies with the physical characteristics of the medium and the zenith angles 
of radiation source and observation. Quantitative calculation of these two components will provide 
important information toward the understanding of the reflective process of solar radiation. 
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Snowpack is a complex and highly variable dielectric medium (References 6 and 7). It is continu- 
ously undergoing metaniorphism due primarily to Interaction with the atmosphere, Direct deposition 
of water vapor on the snow surface will lead to the formation of large (about 0.5 cm) faceted grains 
called surface hoar. Under constant below-freezing temperatures, the grains arc generally well- 
rounded and tljeir sizes vary from a fraction of a millimeter to about one millimeter. The effect of 
a temperature gradient will produce oblong and faceted grains whose sizes vary from about a milli- 
meter to a few millimeters. Tlie melt-freeze process leads to irregular globular masses of polycrystal- 
linc grains and the formation of surface crust, A model of snow will be a medium of randomly 
placed ice grains suspended in air and connected by intergranular bonds, 

Tire “snow surface” is very different from the surface of other media (Reference 8). It could be 
considered as the diffuse interface between two media with different refractive indices. Large scale 
surface rougliness may develop due to wind action. Ripples and wavelike structures are commonly 
observed in prairie snowpacks. Surface slope variation will cause nonuniform solar heating and con- 
sequently nonuniform metamorphism. Local melting and fusion are expected to cause horizontal 
inhoinogenicty. The commonly observed glittering of snow on sunny days is vivid evidence of the 
nhrror-Iikc reflection by microareas of the surface. 

The reflectance of snow shows a strong spectral dependence; moreover, it changes with solar eleva- 
tion and cloudiness (References 9 and 10). Tlie reflectance is high (0.9) for visible radiation and low 
(0.3) for near-infrared radiation. Although the shape of the reflectance curve does not change with 
the advancing of snow metamorphism, the magnitude for different portions of the curve does vary. 

At high solar elevations snow reflects as a diffuse reflector, but at low elevations it showed specular 
characteristics (Reference 9). For a completely overcast condition, snow surface behaves like a 
diffuse reflector (References 10 and 1 1). 

Recent developments in remote sensing via visible and near-infrared radiometers provide great oppor- 
tunity to gather information over snow covered areas. In fact, a data base of snow observations 
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collected by different satellites already exists. These observations have been made at different 
solar zenith angles and view angles as well as different spectral ranges. Understanding the basic 
reflective properties of snow is nccderl in order to analyze both the currently available and forth- 
coming observations. These analyses indicate the calculation of hemispheric reflectance (j.e,, radia- 
tion reflected in all directions) from a limited number of angular observations and calculations of the 
total solar energy reflected from a limited number of spectral band observations (Reference 9, 12, 

13 and 14). 

Tlie theoretical understanding of solar reflectance of snow should be based upon a model which can 
explain the various aspects of its spectral and angular reflection characteristics. In this respect, the 
existing models are not completely satisfactory. Existing models can be classified into two categories. 
One set of models attempts to explain the spectral reflection characteristics, while the other set ex- 
plains the angular variations (References 8, 10, 15 and 16). A formulation therefore is needed to 
provide a unified description of both the spectral and angular variations of reflectance. Tliis paper 
attempts to provide a unified theory for the reflection characteristics of snowpack. Numerical results 
derived from the model are presented and compared with observations, 

II. Formulation of the Reflectance 

The reflective processes of a dielectric medium depend upon surface and subsurface scattering and 
absorption characteristics. The theory of the angular variation of reflected radiation developed in 
Reference 8 is based on the theories developed to explain the reflection characteristics of unpolished 
dielectric surfaces (see also Reference 17). The surface is modeled as a collection of randomly 
oriented facets, each of which reflects like a small mirror, Tlie angular distribution of the incident 
collimated radiation, as reflected by this collection of facets, depends directly on how the mirrors 
are oriented with respect to the mean smooth surface. The understanding of spectral reflectance is 
based upon surface and subsurface multiple scattering by ice grains (References 15 and 16). 
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Scattering and absorption by ice grains depend upon their size and shape as well as the refractive 
index of ice. Due to spectral variation of the retractive index, the scattering and absorption coeffi- 
cients (and hence tlje reflectance itself) vary with wavelength. Application of the radiative transfer 
theory to this multiple scattering problem led to a spectral reflectance which agreed well with ob- 
servations (Reference 16). These two different processes will be discussed in the following sections. 

(1) Surface Reflection 

The intensity of solar radiation on the snow surface can be written as: 

F F 

Ii(O,0) = - g 6 (m - + (1 - g) — (2-1) 

fi 7T 

where 

Fp = total spectral flux incident on the snow surface 
g = ratio of direct and total spectral fluxes incident on the surface 
Mo = cos Oq 
M = cos 0 

(j)= azimuth angle of direct radiation. 

Tlie first term on the right side of Equation (2-1) is the direct radiation and the second term is the 
diffuse radiation. Subscripts for wavelength dependence will be omitted to simplify the notation. 

At the snow surface, part of the direct radiation will be reflected back into the atmosphere. If 
fr(0, 0; O', 0') denotes the bidirectional reflectance distribution function of the surface, then the 
intensity Ij- and the flux F,. of the reflected radiation may be obtained from the following 
equations: 


Ir(0'.0')=^ 


^‘dMM fr(0, 0; O', «A')I(i(0, 0) 

(2-2) 



/*! , , , 




f dM'M'lr(<^',0) 
0 

(2-3) 
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where 


fi' ■ cos O' 


0) “ 7T g 5(/i - Ii^)d(il>) 

Note that following previous formulations (References 15 and 16), diffuse radiation is not considered 
in E<iuation (2-2), For diffuse radiation, the surface reflection is less than 4 percent. Also, it will be 
showji that the surface reflection of direct radiation becomes important only at low elevation angles. 
The bidirectional reflectance distribution function of a typical surface can be illustrated by two well- 
known results. For reflecting surfaces obeying Utmbert’s law (a diffuse reflector) and geometrical 
optics (a specular reflector), the bidirectional reflectance distribution functions are, respectively, 


L(0, ^ ; O', 0') = — (diffuse reflector) 

* 7T 

piO) 


f^(d, <l>; o', <ji') = — • S(p - pi')5(ct>' - 0 + 7t) (specular reflector) 

where a is tiic albedo of the rcllecting surface, and p(0) is the Fresnel reflectivity given by 

rnr--L sin- 0 - cos 0 ^ ^ cosO - \/e- sin^ Q 

^ 2 LVv/*^” sin- O + coso) \ e cos d + - sin^ T/ J 


(2-4) 

(2-5) 


( 2 - 6 ) 


where e is the dielectric constant for the reflecting element. The surface of snow has been modeled 
as a rough surface consisting of randomly oriented tneeted grains (Reference 8). The rough- 
ness scale of the surface can be assumed to be large for visible and near-infrared radiation. Uneven- 
ness of the surface, such as snow on the ridges of a hill or ripples on a prairie snowpack and local 
ice crust, have not previously been considered. These factors are expected to influence the observa- 
tions taken from large snowpack areas. The effect of surface topography on tlie reflected radiation 
can be studied througli the bidirectional reflectance distribution function. To use the rough surface 
scattering theory results developed for radar cross-section analysis (Reference 1 8), the bidirectional 
reflectance distribution function defined in Equation (2-2) is related to the differential cross-section 
per unit area, Y(0, (j>; O', (j>'), generally obtained in radar analysis as: 


^ ' '■ 


( 2 - 6 ) 


where 


U0,r,0\4>')' 


1 

47T cos o' cos 0 


<k') 


1 

cos 0 cos o' 


SIEoP 


(2-7) 


Eg “ scattered amplitude at distance R fr'^m scattering center 
S » surface area of scattering 
Eq » amplitude of incident radiation. 

If there is no preferential orientation of the facets and their distribution is Gaussian (forming the 
so-called Gaussian isotropic surface), then from Reference 19 the bidirectional reflectance distri- 
bution function is obtained as: 


scc'^aexp 






tan^a 
2s^ - 


SiO,0')pW 


Stts^ cos 0 cos o' 


( 2 - 8 ) 


where 


= variance of surface slope 
[sin^ 0 -f sin^ O' ~ 2 sin 0 sin O' cos (j)'] 


tan oc - 


(cos 0 + cos O') 


cos i/' = —^(1 -f cos d cos d' - sin 0 sin O' cos 
v2 

p{\l/) ~ Fresnal’s Reflectivity for angle \{/. 

S(p, o>) = i 

Cq + Cj -f I 


o /2s2f 


tan 0 _ g^fc 


\ 

2cj = tan O' - erfc 

erfc = complementary error function 
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This reflcctanw distribution function resembles one previously used for snow surfaces 
(Reference 8). 

The surface parameter s^ may not be considered an adjustable paramettf lixperimental determina* 
tion of tin's parameter is discussed in Reference 20. From liquation (2-2) and (2-3) the intemsity 
and tlu.x of direct solar radiation rellected by the surface ice grains can be written as: 

V?') * boSM^o* <A') C2-9) 

F,.(do)*Fogf(0y) (2-10) 

where the directional surface reflectance ROq) is defined as: 

d0'jT d/i'/a' f,.(dQ,0;d', 0'). (2-11) 

This surface reflection is primarily responsible for angular vaiiation of radiation reficctcd by the 
snowpack , 

(2) Volume Scattering Ilffect 

Incident solar radiation not reflected by the surface will penetrate through the snowpack, getting 
weaker due to scattering and absorption by the ice grains. As a result of multiple scattering, a frac- 
tion of the incident flux will be reflected back into the atmosphere. The total incident flux for 
which the multiple scattering efiect will be studied is given by (from Equations (2-1) and (2-10)) 

F==Fog[I-f(Oo)]+(l-g)Fo 

= Fo[l~gf(d^,)] (2-12) 

This effect may be examined by studying the solution of the radiative transfer equation (Refer- 
ence 21). This is the energy conservation relation for radiation propagating through a scattering and 
absorbing medium. The spectral absorption coefficient of ice, the grain size, shape, density and 
thickness of the snowpack enter into this equation through parameters called the single scattering 
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albedo w, the extinction coefficient % the optical thickness n and ihe phase function p. The inten- 
sity of scattered radiation within fsc snowpack may be obtained from ine solution of the following 
equationj 

M * “Kr, m) y'j p(ti, /i') Ur, iA')dn' P(/<. IJ^q) (2-13) 

where 

I * intensity of radiation 

f 7(x')dx' 

0 

X » depth within snow from the surface. 

The single scattering albedo (co) is defined as the ratio of scattering and extinction cross-sections of 
the ice grains. For this model the Mic scattering theory will be used to calculate the single scattering 
albedo (Reference 23). A simple formula which reproduces the Mic theory results for large spheres 
with reasonable accuracy (References 22 and 23) is 

w®| + ^exp(-ckr) 

where 

c= 1.67 

k = spectral absorption coefficient of ice (References 23 and 24) 
r = radius of ice grains . 

Although the Mie scattering calculation is based on spherical particles, it has been suggested that 
Equation (2-14) may be u.scd for nonspherlcal particles by choosing a different value for the 
constant c (Reference 22) or by using the radius value of a circular area equivalent to the shadow 
area of the particles (Reference 25). 

The extinction coefficient (T) is the product of the extinction cross-section ( oq ) and the number 
density (N) of ice grains. For large spherical particles one obtains: 


(M5J 


a„ ■ 2itr 


4irP \Pi/ 


i2Ab) 


where r is the radius, and and p\ are respectively tlie densities of snow and ice. TIuis the extinc- 
tion coefficient is jjiven by 


m ■ 


(2-17) 


The scattering phase function. y>(p,p‘), which is the ratio of differential and total scattering cross- 
sections can also be obtained from the Mie scattering theory. Since large i)articles scatter strongly in 
the forward direction, the j<uase function is strongly asymmetric (Reference 15). Further details 
about this function will be discussed later. 


The boundary conditions for the solution of this equation are: 

At the snow-air interface 

i+(r=0,M) = 0 (2-18) 

and at the bottom interface, assuming soil can be treated as a Lambertian reflector, tlie boundary 
condition is 

a^F 

L(r = Tq, p) = I+(r = /u) +~ exp (-Tq/Pq) (2-19) 

where 

Hg = albedo of soil surface 

Tq ™ optical thickness of the snowpack 

and, 4 and L are respectively the intensities going away and towards tlie surface. 

Accurate solution of the radiative transfer equation for strongly asymmetric phase functions is nu- 
merically difficult and time consuming. However, this asymmetric nature of the phase function has 
been used to develop approximations which are reasonably accurate (References 2b and 27). Tlie 
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idea is to approximate the forward scattering peak by a delta function, Once the forward scattering 
peak is subtracted from the phase function, the remaining scattering function is amenable to efficient 
and accurate solutions. By writing the phase function as 


p(M, ix') = 2d2 6(/u - /r') + (1 - fx') (2-20) 


(where is the fraction of radiation scattered in the forward direction, and p*(fx, fx') is the peak 
subtracted normalized scattering function) one can show that the form of the radiative transfer 
equation remains invariant. The transformed radiative transfer equation contains, in place of 
parameters r, to, and p(/a, fx'), the new parameters r*, w'*‘, and p*(fx, p'), where 


T* = (1 - COd^)T’ 


( 2 - 21 ) 


(I-d^)co 
(1- d^w) 


(2-22) 


The value of d^ for large semitransparent ice spheres is approximately 0.76 in the visible region, and 
0.92 in the near-infrared region. From above discussion, these new parameters are to be used in 
Equations (2-13), (2-1 8) and (2-19) to study the volumetric scattering effect. 


(3) Spectral Reflectance and Albedo 

If the spectral reflectance, A\, is defined as the ratio of reflected and incident fluxes of mono- 
chromatic radiation, then one obtains 


Ax(0o) = 


Fogf(0o) + 27 t 


X 


H)fx dp 


= F(0o)g+[l-gf(Oo)l 


(2-23) 

(2-24) 


where I_(0, p) is the intensity of volu metrically scattered emerging radiation, and 

27t /* l_(0,p)pdp 

Jo 


Rh=' 


Fo[l-gf(0o)) 


(2-25) 


is the diffuse reflectance of snow. The first term on the right-hand side of Equation (2-24) is due to 
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I surface reflection, and the second term is due to multiple scattering. Changes in surface texture and 

topography will affect the spectral reflectance for clear sky conditions, Since the ratio of direct 
and total radiation (g) depends upon the wavelength, the effect of surface texture and topography is 
not expected to be same for all wavelengths. From the spectral reflectance, the albedo can be calcu- 
lated as; 


[ 

' A(0o) 


/"axWo)Po <1^ 

FodX 



(2*26) 


III, Illustrative Numerical Results 

The directional surface reflectances calculated for the bidirectional reflectance distribution function 
equation (2-8) are given in Tabic 1 for the surface parameter s^ = 0.01, The calculation was per- 
formed using 1° azimuth grid with the Simpson’s rule, and the 64 point Gauss quadrature for the 
elevation angle. 


Table 1 : Directional Surface Reflectance 


Solar 

Elevation 

(Degrees) 

5.0 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 

40.0 

Surface 

Reflectance 

0.3403 

0,2422 

0.1739 

0,1273 

0.0873 

0.0622 

0,0458 

0.0355 


From the solution of the radiative transfer equation, the subsurface reflectance for a homogeneous 
snowpack was calculated using the following equation (Reference 28, chapter 8, and Reference 27): 

2 - 3qMo + (1 - OJ*) 3a*Po(2)u - q) 

Rd= (3-1) 

2 (l+qXl-e^Mo^) 

where 

a* = d/(l +d) 

5 = [3(1- cj*)(l - oj*a*)]i/2 
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q»26/f3(l -w*a*)l 


In order to calculate the spectral reflectance of snow, it is necessary to obtain the incident flux of 
the solar radiation. For clear sky conditions, the solar radiation incident on a snowpack can be 
decomposed into two parts: the direct radiation and atmospherically scattered radiation. Tlie 
direct radiation incident on the surface differs significantly from the radiation incident on top of the 
atmosphere. Atmospheric constituents such as ozone, oxygen, carbon dioxide, and water vapor are 
spectrally selective absorbers of solar radiation. The relative magnitude and angular variation of at- 
mospherically scattered radiation varies with solar elevation, atmospheric conditions and the surface 
reflectivity, 

The spectral flux outside the atmosphere (Fqj^) is multiplied by parameterized attenuation coeffi- 
cients of ozone (ro^), water vapor mixed gases (Tg), aerosol and the Rayleigh scatter- 

ing (^■r) given in Reference 29 to obtain the direct incident flux on snow as 

Px ^ = ^0X^03 ’’HjO Va^'R “ ( 3 - 2 ) 

where <x is the solar elevation angle. The incident atmospherically scattered radiation is calculated 
assuming snow albedo to be about 0.75 by the method suggested in Reference 30 from the equation 

r px^ 1 

Pf = £ [PoX’-Oj’-HjO’-g - 37 ; J « (3-3) 

where 

^ = 0.63 for a <7.5° 

= 0.602 + 0.0037a for a > 7.5° 

From Equations (3-2) and (3-3) one can calculate the total flux Fq, and the ratio of the direct and 
total fluxes g as: 

Po + 

g = F^)/Fo 
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An equidistance grid of 45 points from 0.30 to 2.50 micron radiation is chosen to represent tlie inci- 
dent flUK. Using this representation, the spectrally integrated values arc obtained by the Simpson 
rule method. The atmospheric parameters for a dry clear winter condition (References 29 and 30) 
arc chosen as follows: 

Ozone Path Length = 0,28 cm NTP 
Water Vapor Path Length - 0.9 cm 
Angstrom Aerosol coefficient = 0.08. 

The calculated spectral incident flux for two different solar elevations are shown in Figure 1. Ab- 
sorption of ultraviolet radiation by ozone, and near-infrared radiation by water vapor and carbon 
dioxide is evident in the incident solar flux. At low elevations, the redness of sky can be seen through 
the broadening of incident visible flux due to atmospheric scattering. 

The spectral reflectance of snow calculated from Equation (2-24) is shown in Figures 2 and 3, The 
chosen radii (r) are typical snow crystal sizes for different jiages of equi-temperature metamorphism. 
Tlie overall shape of the reflectance curve does not seem to change with the crystal size or the solar eleva- 
tion, but relative magnitude of various parts of the curve does vary. The spectral reflectance in the 
near-infrared regions shows maximum sensitivity to the crystal size and the solar elevation. A de- 
crease in the solar elevation or the crystal size appears to have a similar effect on the spectral rellcc- 
tance. Calculated reflectance is in qualitative agreement with the obseivation (Reference 31). The 
spectrally integrated directional reflectance (the albedo) obtained from Equation (2-26) is shown in 
Figure 4. Quantitative agreement between theory and averaged antarctic observations by Liljeqiiist 
(see Reference 32) is good, but should be considered tentative because a direct comparison with the 
incident solar flux and the spectral reflectance cannot be made. Furthermore, corresponding to 
these observations the snow parameters are not exactly known and the chosen values of the atmo- 
spheric parameters may not be representative of subpolar regions. Note, jiowever, for clear sky con- 
ditions the albedo values are more sensitive to snow parameters than the atmospheric parameters. 
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Figure 3. Spectral Reflectance of Snow for Different Solar Elevation (Grain Radius = 0. 1 mm) 


16 









In this study ‘microscopic’ surface texture (shape and orientation of surface crystais) is modeled 
tlirougli tlie parameter s in liquation (2-8). For faceted grains witir more or less Iiorizontal surfaces 
the vaiue of this parameter is smaller than that for randomly oriented grains. As the value of S de- 
creases, the magnitude of surface reilectance increases i.c., the specular surface reflection increases. 
The slope of the albedo - elevation curve at low elevation angles (for elevation angles less than about 
35 degrees) depends primarily upon tlie specular surface reflection and the grain size. The slope will 
increase with tiie grain size, and higlier specular reflection. At liigli elevation angles, the importance 
of surface reilection is much less tlian the volumetric renection. Consequently, for these elevations, 
tlie albedo depends mainly on the grain size. The slope at low elevation angles for an old metamorpiv 
osed snow with surface crust is expected to be higlier than that for a wind blown snow with rounded 
grains. Thus observed differences in the elevation dependence of the albedo (References 9, 12, 13, 
and 33) may liave been due to differences in the surface texture and the grain size. A quantitative 
comparison with tliese observations would be argumentative because pertinent snow parameters are 
not known, and a comparison with the spectral reflectances cannot be made. A more complete set 
of observations is desirable. 

IV. Conclusion 

A model for solar reflectance of a snowpack is formulated based upon surface reflection and 
multiple scattering by the ice grains. The model calculations show that the surface scattering plays 
an important part in tlie angular variation of reflectance particularly for low solar elevation angles. 

At high solar elevation, the albedo primarily depends upon tlie grain size. Model calculations are in 
qualitative agreement with the observed spectral and integrated rcflectaiKCs, Further study of this 
problem should be directed towards a more exact solution of tlie radiative transfer equation in order 
to study the angular dependence of multiple scattered radiation. For a better understanding of the 
reflective process, a more complete set of observations is needed. 
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